Super- Arrhenius slowing down in the metastable liquid of hard spheres 
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It is demonstrated that a straightforward extension of the Arrhenius law accurately describes dif- 
fusion in the thermodynamically stable liquid of hard spheres. A sharp negative deviation from 
this behaviour is observed as the liquid is compressed beyond its stability limit, which can be com- 
pared with super- Arrhenius slowing down in supercooled liquids. It is suggested that the observed 
dynamical transition is caused by extra entropy barriers arising in the free-energy landscape. This 
interpretation is corroborated by the observation of characteristic anomalies in the geometry of a 
particle trajectory. 

PACS numbers: 66.10.Cb, 64.70.Pf 



Computer simulations of liquids using hard spheres 
(HS) remain a key source of information about the most 
fundamental aspects of the liquid state. The main in- 
tellectual attraction of these models is in their purely 
geometric nature: all the observable properties of a HS 
liquid can be reduced to a single thermodynamic quan- 
tity - entropy a measure of finite-range structural corre- 
lations §,§. Thus, mapping a phenomenon exhibited by 
liquids with continuous interactions on a HS model yields 
its interpretation in terms of the statistical geometry of 
sphere packing. 

An important aspect of liquid behaviour addressed us- 
ing HS models is supercooled liquid dynamics j^, [|, ^) . 
This notion refers to the complex of dynamical anomalies 
observed in a liquid which remains in a metastable equi- 
librium below its melting point. For the one-component 
HS system, the stable liquid domain is bounded by the 
critical value of the packing fraction rj c = 0.493 || 
(rj = it pa 3 / 6, where p is the density and a is the HS 
diameter). Having been compressed beyond r\ = 0.54, 
the HS liquid inevitably crystallizes [||. 

The most prominent indicator of the supercooled liq- 
uid regime, and, in fact, its defining feature ||, is 
the characteristic behaviour of the relaxation dynamics 
which slows down with cooling progressively faster than 
it can be inferred from the Arrhenius prescription univer- 
sally describing the temperature variation of the trans- 
port coefficients in stable liquids. 

In contrast to atomistic models with continuous inter- 
action, the HS systems lack the energy scale. Therefore, 
in order to compare the non- Arrhenius liquid dynamics 
as observed in the former with the respective behaviour 
of the metastable HS liquid, we have to resolve the funda- 
mental problem of finding a description of the stable HS 
liquid dynamics that would be an adequate analog of the 
Arrhenius law. In this capacity, Hildebrant-Batchinski 
relation || connecting the diffusion coefficient with the 
free volume was suggested [|l0|. The HS liquid diffusion 
Jl0| follows this relation in the stable domain and de- 
viates from it for r\ > r\ c W%- However, this deviation 



is positive and, therefore, cannot be regarded as an ap- 
propriate analog of the super- Arrhenius slowing down in 
supercooled liquids. 

In this Letter, it is demonstrated that an earlier sug- 
gested relation between the diffusion coefficient and the 
thermodynamic entropy [jl2| can be interpreted as a 
straightforward extension of the Arrhenius law. Hav- 
ing been tested by the molecular dynamics, it is found 
to accurately describe the diffusion in the stable one- 
component HS liquid. Moreover, a significant slowing 
down as compared with this relation is observed in the 
metastable HS liquid domain, a direct analogue of the 
super- Arrhenius slowing down in supercooled liquids. 
This transition in the HS dynamics is discussed in terms 
of the free energy landscape. 

The Arrhenius law asserts that the diffusion coefficient, 
D, scales with the temperature T as 



D = D e- A/kE 



(1) 



where A, the so-called activation energy, is interpreted as 
the average height of the free-energy barrier, per particle, 
that has to be crossed in order to perform an elementary 
step in the process of diffusion. It can be presented as 
A = E a — ksTSa, where E a and S a represent, respec- 
tively, the heights of the energy barrier and the entropy 
barrier involved. For the HS liquid, A/ksT = —S a ; e Sa 
can be understood as the inverse average number of at- 
tempts per a successful diffusive transition of a particle. 
Two possible types of the free-energy landscape can be 
discussed in the context: 

(i) The available configurations are abundantly con- 
nected; the connectivity is facilitated by independent mo- 
tions of individual particles, and each particle has an un- 
restricted immediate access to the entire space volume 
assigned to it by the structural constraints. Then the 
inverse average number of attempts needed for a particle 
to move from its current position scales as e _j4 / fcsT = e s 
where s is the excess entropy, per particle, - the difference 
between the system's entropy and that of the perfect gas 
at the same thermodynamic conditions. 
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(ii) The configuration-space connectivity is re- 
stricted by additional constraints, complementary to 
the ensemble-averaged structural constraints dominating 
type (i) landscape. In this case, the diffusive motions of 
an individual particle are coupled, in a presumably hi- 
erarchical manner, with the respective motions of other 
particles within a certain range [ ft3| . These dynamical 
constraints give rise to extra high entropy barriers which 
divide the configuration space into a set of components 
jl4j]. These barriers dominate the relaxation dynamics 
and, in this way, control the diffusion rate. In terms of 
Eq. (1) this means that e ~ A / kBT < e s It has to be em- 
phasized that this type of free-energy landscape exists on 
a limited time-scale and, therefore, is quite distinct from 
the equilibrium distribution. However, there exists a con- 
jecture fl5|j expressing the cooperativity range in terms 
of equilibrium entropy; this conjecture is not discussed 
here. 

The described landscape regimes can be understood by 
keeping in mind an important topological peculiarity of 
the multidimensional configuration space where each 
point is adjacent to almost all other points. This can be 
illustrated by the "small world" topology of abundantly 
connected multidimensional networks fL7j , On the other 
hand, if an entropy barrier characteristic of type (ii) land- 
scape arises, it will be immediately close to almost all 
points of the configuration space with a profound effect 
on the general rate of dynamics. 

The central postulate we adopt here is that the free- 
energy landscape that controls the relaxation dynamics 
in the stable HS liquid is of type (i). Having assumed this, 
Arrhenius law (1) can be transformed into the following 
simple form relating s and D: 



D = D Q e s 



(2) 



If the time-scale is measured in terms of the Enskog col- 
lision frequency T |H, this relation becomes equivalent 
to the earlier suggested scaling law for atomic diffusion 
|jj"2"f which was successfully tested on a number of simple 
liquids, including the HS liquid using a two-body approx- 
imation for s pi. In the full s version, the scaling law was 
tested on some metallic liquids but for the HS liq- 
uid, such tests have not been done so far. In this context, 
a conjecture should be mentioned where D oc e Bs with 

In this study, relation (2) was tested using a molecular- 
dynamics simulation of a HS system comprising 6912 
identical particles. This simulation was necessitated by 
the fact that, although the HS diffusion has been ex- 
plored in a number of simulations, a considerable uncer- 
tainty remains concerning its size-dependence at different 
densities §, 

The assumption that the dynamics in a liquid is con- 
trolled by the excess entropy, a measure of the struc- 
tural correlations, implies that the length should be mea- 
sured in terms of the characteristic length of the structure 
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FIG. 1: The structure factor and the radial distribution 
function of the HS liquid. Solid line, 77 = 0.49; dash-dotted 
line, 77 = 0.38. The dashed line, in both plots, the Lennard- 
Jones liquid simulated at p — 0.85 and T = 0.7. 



which manifests itself in the position of the main peak of 
the structure factor. Fig. 1 shows that the latter scales 
with the density p as p 1 ^ 3 ; therefore, p -1 / 3 is a conve- 
nient unit of length. Fig. 1 also shows that this unit 
of length is consistent with the equilibrium separation 
of the nearest neighbours in the Lennard- Jones liquid of 
which the structure factor can be reproduced by that of 
the HS liquid. The use of a as the unit of length that was 
adopted in the earlier test of the scaling law Jl2| on the 
HS liquid is apparently inappropriate at lower densities. 

The HS excess entropy was calculated here by the ther- 
modynamic integration using Boublik-Nezbeda approxi- 
mation p3| for the equation of state. A minor correction 
was introduced for 77 > 0.47 ||. 

The simulation results are shown in Fig. 2. Each value 
of D was obtained from the mean-square displacement 
averaged over 10 4 collisions per particle. These results 
are compared with another set of the HS diffusion data 
which represent an extrapolation to the infinite system 
size|2l], g^]. It is clear that with Dq — 0.079, relation (2) 
accurately describes the diffusion in the one-component 
HS system within the range 0.36 < 77 < 0.49 which covers 
the entire stable liquid domain. 

Do corresponds to the diffusion rate in a hypothetical 
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FIG. 2: Diffusion coefficient in the HS liquid as a function of 
the packing fraction r\. Dots, the present molecular dynamics 
simulation; dotted line, the extrapolated limit for the infinite 
system size j2l], chain-dashed line: Eq. (2)XS with Do = 
0.078. Dashed line: the Enskog approximation pM. 



situation where the discussed liquid diffusion mechanism 
is still valid but s = 0. The observed value of Dq can 
be rationalized using the following arguments based on 
a simple toy model in the spirit of the lattice gas mod- 
els. Consider a tessellation of the space into equal-size 
cells. Each cell contains one stationary particle which 
stays within the cell until its status is changed; besides, 
there is a number of mobile particles distributed in space. 
At each time-step, a mobile particle collides with the sta- 
tionary particle that is in the same cell. As a result, the 
former becomes stationary, whereas the latter becomes 
mobile and traverse the cell boundary into a randomly 
chosen adjacent cell (if several mobile particles simulta- 
neously occur in one cell, an arbitrary collision sequence 
is assumed). Thus, each particle performs a random walk 
with the step size rj p -1 / 3 . Assuming that the number of 
mobile particles is small, and keeping in mind that there 
is one particle jump per two collisions, the diffusion co- 
efficient can be estimated as D w Tp^ 2/>3 /12, which is in 
a good agreement with the above value of Do. 

The results presented in fig. 2 demonstrate that in 
the metastable liquid domain r\ > ?/ c , the HS diffusion 
exhibits a rapidly increasing slowing down as compared 
with the Arrhcnius-like behaviour conjectured by relation 
(2). Apparently, this behaviour can be regarded as a 
direct analog of the super- Arrhenius slowing down in the 
supercooled liquids with continuous interactions. Using 
the above arguments, the onset of this new dynamical 
regime can be interpreted as a clear signal that the height 
of entropy barriers exceeds the ensemble-averaged excess 
entropy, and the free-energy landscape that controls the 
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FIG. 3: The ratio of the mean square maximum displacement 
(Z 2 {t)) to the mean square displacement (R 2 (t)). Solid line, 
r\ — 0.48; dashed line, rj = 0.52; chain-dashed line, r\ = 0.53. 
Dotted line indicates the random walk limit (Z 2 (t)) / (R 2 (t)} = 
1.4024... 0. 

HS liquid dynamics transforms from type (i) to type (ii) . 

The long-time configuration-space constraints con- 
jectured to dominate the HS liquid dynamics in the 
metastable domain can be detected from their impact 
on a particle trajectory. The latter represents a 3D real 
space projection of the multidimensional trajectory of the 
system in its configurational space. The relaxation dy- 
namics that unfolds in type (i) landscape represents a 
random walk in the configurational space, which must 
also be true for a single particle trajectory. The long- 
time constraints caused by the development of extra bar- 
riers characteristic of the type (ii) landscape give rise to 
distinct anomalies in the geometry of a particle trajec- 
tory. A quantity that can be conveniently employed as an 
indicator of a possible deviation from the random walk 
geometry in a particle trajectory is the maximum abso- 
lute value of the displacement from the original position 
within the time interval (0,t): 

Z(t) = max{R(r), < r < t} (3) 

The ratio of its second moment {Z 2 (t)), the mean square 
maximum displacement, to the mean square displace- 
ment (R 2 (t)) for a random walk in 3D space is equal 
to 1.4024... p4| . The evolution of this ratio as a func- 
tion of (R 2 (t)} for various values of r\ is presented in 
Fig. 3. For 77 = 0.48, in the stable liquid domain, the 
random-walk limit is attained after a short initial pe- 
riod of ballistic behaviour, as soon as a particle leaves 
the cage of its nearest neighbours. The geometry of a 
particle trajectory changes significantly as the HS liquid 
is compressed beyond its stability limit: a new regime of 
diffusion develops where (Z 2 (f)) / (R 2 (t)) exceeds the ran- 
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dom walk value. This indicates that a particle trajectory 
has a higher chance than random to return to an earlier 
covered region. The extent of this apparent confinement 
effect in a particle trajectory amounts to the diffusive dis- 
placement of several particle diameters, and it increases 
rapidly, as well as the magnitude of the effect, with the 
increase of rj. 

According to the arguments presented above, the 
height of an entropy barrier can be regarded as a mea- 
sure of the real-space extent of the local configurational 
transformation that the system has to perform to cross 
the barrier. Extra high barriers in type (ii) landscape 
correspond to extended correlations, both positional and 
dynamical. The rapid increase in the (time-limited) cor- 
relation length in the metastable HS liquid that has been 
detected here from the geometry of a particle trajectory 
represents a generic feature of the supercooled liquid be- 
haviour p5| . Ergodicity restoring relaxation whereby the 
extra high entropy barriers are created and destroyed is 
facilitated by highly collective activated hopping. This 
dynamics was discerned in the metastable HS liquid us- 
ing dynamical density-functional theory analysis p(| . 

The profound similarity between the pattern of slow- 
ing down in the metastable HS liquid and that in con- 
ventional supercooled liquids indicates the entropic na- 
ture of this phenomenon p^] , and suggests a unifying 
scenario for its development in terms of the free-energy 
landscape transformation. The change in the packing 
geometry under cooling/compression beyond the liquid 
phase stability limit, presumably caused by reduction of 
the free volume, leads to formation of long-lived extra 
high entropy barriers and the loss of the short time-scale 
configuration-space connectivity. This results in chang- 
ing the relation between between s and D which can be 
detected as a deviation from the scaling law (2). 

This study was supported by the Swedish Research 
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